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Abstract 

With Monte Carlo simulations, we investigate the relaxation dynamics with a domain interface 
in a magnetic system at zero temperature, taking two-dimensional driven random field Ising model 
with quenched disorder as an example. The dynamic scaling behavior is carefully analyzed, and a 
dynamic roughening process is observed at the depinning transition. For comparison, additional 
simulation without disorder has been performed as its background. The effect of the overhangs is 
discussed, and the growing interface exhibits intrinsic anomalous scaling and spatial multiscaling 
with Q > 1 and Cloc < lj different from the work described by QEW. 
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I. INTRODUCTION 



In the past years much progress has been achieved in the study of dynamic processes far 
from equilibrium. For example, the universal dynamic scaling form in critical dynamics has 
been explored up to the macroscopic short-time regime jl-4|, when the system is still far 
from equilibrium. Although the spatial correlation length is still short in the beginning of 
the time evolution, the dynamic scaling form is induced by the divergent correlating time 
around a continuous phase transition. Based on the short-time dynamic scaling form, new 
methods for the determination of both dynamic and static critical exponents as well as the 
critical temperature have been developed J2-5]. Since the measurements are carried out in 
the short-time regime, one does not suffer from critical slowing down. Recent progress in 
the short-time critical dynamics includes, for example, theoretical calculations and numerical 



simulations of the XY models anc 
disorder jz ioj], ageing phenomena 



Josephson junction arrays [6], magnets with quenched 
111 ], weak first-order phase transitions , and various 



applications and developments 13]. 

Meanwhile, the physics of elastic systems in disordered media also has been the focus of 
intense theoretical and experimental studies in the different areas, such as charge density 



waves 



Q, 



materials, contact lines 
[2 



20 



vortex lattices 15 L domain walls in magnetic 16|, Il7 | o r ferroelectric 



2l| . fluid invasion of porous media 



18 



19] 



22], dislocations moving 



[ [24]], and crack propagation [25j. With increasing driving force F, the driven interface 
in quenched disordered system displays a transition from a pinned interface to a moving 
interface. This so-called pinning-depinning phase transition at zero temperature is viewed 



as a critical phenomenon 



26N35| . and its ordered parameter is the interface velocity v. 



With the assumption that the interface shows the properties of an elastic membrane, this 
growth process can be described by the driven Edwards- Wilkinson equation with quenched 



32- 



or the driven quenched Kardar- 



disordered (QEW) in the isotropic universality class 
Parisi-Zhang equation (QKPZ) 39M42| in the anisotropic universality class. With these 
equations, all the critical exponents can be derive d by simulation and theory analysis called 
functional renormalization group (FRG) 29|, |36|, |43| . And the roughness exponent ( > 1 



obtained from the equation QEW in 1 + 1 dimensions means that the elastic string will 



inevitably break as the length of the string is increased 



34 



23, 



38 



Recently, these kinetic roughening phenomenons have brought about interest on the su- 
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perrough growth of the interface or surface with ( > 1, not only in simulations 44H46j . 



but also in experiment 



47H49I] . Three scaling behaviors have been observed experimentally, 



namely, (i) FV standard scaling for ( = (i oc = ( s < 1; (ii) intrinsic anomalous scaling for 
C 7^ Cloc — Cs < 1; an d (hi) superrough scaling for ( = ( s > Q oc = 1. Here roughness 
exponent (, local roughness exponent Q oc and spectral roughness exponent ( s can be mea- 
sured in different ways, such as globe width function W(L,t), local width function W(l,t), 
height correlation function C(r,t), local interfacial orientation W n (l,t), and power spectral 
density S(k,t). And spatial multiscaling behavior is observed associated with the intrinsic 



anomalous scaling class (ii) 



44 



|47j. 



On the other hand, many recent activities have been devoted to the domain-wall dynamics 
in the magnetic materials, which is an import ant topic in magnetic devices, nano-materials, 
thin film, and semiconductors 



17 



5CH54I] . And external field H or electrical current J 



is applied as the driving force F. At zero temperature, the depinning transition is obtain 
when the homogeneous driving field H reaches its critical value H c , which is relation with 
the disorder. And if a periodic external field H(t) = Ho cos(o;t) is applied and/or a non-zero 



temperature is introducec 



phase transitions 



52 



53 



.t 



re domain wall exhibits different states of motion and dynamic 



55 



561 ] . Most these works concentrate on the stationary state at 



the zero or low temperatures and in response to the external magnetic field H(t). Instead 
of the phenomenological model, such as QEW or QKPZ equation, a well established model 



to investigate the dynamics be 



ravior of the domain interface is the driven random field 



Ising model (RFIM) 



26 



28 



57 



591 ] . The scaling behavior and critical exponents have been 



obtained for various dimensions at zero or low temperature, which is roughly consistent 
with the result of QEW. While as we well known, the roughness growth is rarely to be 
investigated because lots of overhangs occur at the depinning transition and it is hard to 
determine the critical exponents and transition point precisely by the finite-size scaling and 
finite-temperature scaling analysis in the stead state 

Very recently, the experiment on the non-stationary dynamic evolution of the domain 
wall is reported for the creep and pinning effect [60|, though its dynamic behavior is de- 
scribed roughly. What's more, the short-time relaxation of a driven elastic string in a 
two-dimensional pinning landscape is also reported in the reference 33j, in order to get 
the critical exponents of the depinning transition independently. Additional, the short time 
dynamic relaxation of a domain wall has been concerned for magnetic systems at a stan- 
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dard order-disorder phase transition in the Ising model jsl, Q, [13(|. It is described by the 
relaxation dynamics starting from a semi- ordered state, and shares certain common features 
with those around free and disordered surfaces. Since no external magnetic fields are added, 
macroscopically the domain wall does not move, but a kind of roughening phenomena occurs. 

In this paper, we investigate the relaxation dynamics of domain walls at zero temperature 
for the depinning transition, taking the two-dimensions driven RFIM as an example. And we 
aim at understanding the non-stationary properties of this dynamics system and determining 
the static and dynamic exponents as well as the transition point. Then the effect of overhangs 
will be studied by the roughness scaling behavior, which causes the main difference between 
QEW and RFIM at the depinning transition. For comparison, its background in our work is 
also investigated without overhangs. In Sec. II, the model and scaling analysis are described, 
and in Sec. Ill, the numerical results are presented. Sec. IV includes the conclusions. 



II. MODEL AND SCALE ANALYSIS 

The random-field Ising model is defined by the Hamiltonian 

^ = " J E SiSj-H^Si-^hiSi. (1) 

<ij> i i 

Where Si = ±1 is just the Ising spin on the two-dimensional lattice with Lx L. The random 
field hi is taken from a distribution within an interval [—A, A] and hence is bounded. H 
is the homogeneous driving field. The strength of the random field is fixed as A = 1.5J in 
this paper, and A = OJ is also used in comparison as its background. For convenience, we 
set J = 1 in this paper. The main simulation is performed with the lattice of size L = 512 
at zero temperature, and the maximum updating time is = 1024. Additional simulation 
with L = 256 and L = 1024 are also performed to exclude the finite-size effect. The total 
samples for average are 50000, and the statistical error are estimated by dividing the samples 
into two or three groups. If the fluctuation in the time direction is comparable with or larger 
than the statistical error, it will be taken into account. 

The initial semi-ordered state is such a state that spins are positive in the sublattice on 
the left side x < 20 and negative on the right side x > 20. Here we set the x axis in the 
direction perpendicular the initial interface and its left boundary is located at x = 0. Then 
all the system is rotated by an angle 7r/4 in order to make sure that the initially interface 
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of the system is just in the (ll)-direction of the lattice. Therefore there is no pinning of 



the interface without the disorder 



26- 



The antiperiodic conditions are used in the x 
direction. While in the other direction, periodic conditions are used. After preparing that, 
we perform Monte Carlo simulations and the spins are chosen randomly to update. And 
the spins will be flipped only when the total energy decreases after flipped. As time evolves, 
the interface moves due to the driving field and becomes more and more roughening. In 
Fig. [D(a), the real snapshot images of this interface are illustrated in the dynamic relaxation 
process for the depinning transition and its background. The local interface at the transition 
point not only tilts, but also forms local grooves or spikes, quit different from the instance 
of the background. Hence, it is difficult to define the height function h(t) of the interface 
because of visible overhangs. 

In this case, we define the velocity of the domain interface as the time derivation of the 
]. First, the magnetization and its second moment are measured, 



magnetization 



M«(t) 



1 

J2k 



x,y=l 



1,2. 



(2) 



Here S xy (t) is a spin at the time t on the lattice position (x,y), and < ■ ■ ■ > represents the 
statistical average over disorder. For convenience, we also use M(t) = M^(t) to denote the 
magnetization. Then the velocity of the domain interface is shown, 

L d M(t) 



v(t) 



(3) 



2 d t 

Here L/2 is the scale factor. In order to characterize the growth of the domain interface and 
its fluctuation in the x direction, we introduce a height function and its second moment 



1 



.x=l 



k = 1,2. 



(4) 



Here < ■ ■ • > represents not only the statistic average but also the average in the y direction. 
As usual, we also use the notation h(t) = h^\t) and the relation h(t) = M(t) is obtained. 
Then the roughness function of the domain interface is defined, 



uj 2 (t) = hP\t)-h{t)h(t). 



(5) 



After that, the relation between the magnetization fluctuation and interface fluctuation can 
be described in the function F(t), 

hfW(t) - M(t)M(t) 



F{t) 



hW(t)-h(t)h(t) 



(6) 



What's more, a more informative quantities for the interface, called the height correlation 
function C(r,t), is introduced. 

C(r,t) = ([h(y + r,t)-h(y,t)] 2 ). (7) 

Here h(y, t) describes the local interface height and can be obtained by Eq. (j4|) with a fixed 

y- 

At the depinning transition point H = H C (A), the main scale length in the dynamic 
system is the non-equilibrium spatial correlation length For a finite size, the lattice 

size L is an additional length scale. And the system exhibits as the second-order transition 
when temperature is zero. So one can believe that £(t) grows as a power law £(t) ~ t x l z and 
z is so-called dynamic exponent. General scaling arguments lead to the scaling form of the 
order parameter v(t). 

Here ft and v are the static exponents and r = H — H c is defined as the departure from 
the transition point. And the scaling function G (£(£)/ X, r/£(i)) is constant if L — > oo and 
r = 0. Then the scaling form is simplified to 

v(t) ~ r? /vz . (9) 

And the scaling behavior of the derivative d T \nv(t,r)\ T=0 can be obtained according to 
Eq. (H. 

d T \nv(t,T)\ T=0 ~t l ^ z . (10) 

In general, the height function h(t), the roughness function oo 2 (t) and the height correla- 
tion function C(r, t) in Eqs. (j4j), ([5]) and ([7]) do not obey a simple power law behaviors. Since 
the depinning transition is a dynamic transition and its order parameter is v (t) ~ dh{t)/dt. 
And in the background without disorder, the interface is also roughening. So pure roughness 
function and height correlation function should be defined by subtracting the contribution 
from the background. 

Du\t) = co 2 (t) -co 2 (t). (11) 
DC(r,t) = C(r,t)-C b {r,t). (12) 
And the scaling form of Duj 2 (t) and DC (r, t) are shown as in the case of a standard growing 



interface, reported in Refs. [30 



31 



44). 



Du 2 {t) ~ r 2 ^ z . (13) 



m)2(c-C, oc ) r 2C ioc if r «£(t) 
DC(r,t)~4 . (14) 

e(t) 2C if « r 



Here ~ t 1//z and z have been denned before. And ( is the roughness exponent and £, 



z ot- 



is the local one. Since the region with the power law behavior 



30 



'or DC (r, t) vs. r is rather 



3l|. 



narrow, a new scaling form is introduced to fit this relation 

DC(r, oo) = A[tanh(r / 'B)] Kl °* . (15) 

Then, the roughness function oo 2 (t) represents the fluctuation only in the x direction, while 
the fluctuation of the magnetization is measured in the whole lattice. So the scaling behavior 
of the function F(t) is shown, 

F(t) = m 4 ' 1 ~ t 1/z - (is) 

What's more, all the scaling behaviors above hold only after a time scale t mic , which is 
sufficiently long in the microscopic sense, but still short in the macroscopic sense. Hence, a 
power law correction to the scaling is used sometimes in order to fit the data well. 

Y(x) ~x a (l + c/x). (17) 

Here the convergence of Y(t) is to a power law behavior and the parameter a is viewed as 
the critical exponent. 

III. THE RESULT OF NUMERICAL SIMULATION 

Firstly, the depinning transition of the 2D driven RFIM is measured carefully by the 
short time dynamics in Fig. QJb) with different driving filed H and the disorder is fixed 
as A = 1.5. The transition point H c = 1.2933(2) is obtained according to general scaling 
theory, which is more precise than H c = 1.290(5) obtained through the stead state near the 



depinning transition [28|. Then the exponent (3/vz = 0.217(2) is measured from the slope 
of the curve at the transition point, according to Eq. ()9]). And the finite size effect is also 
checked with the different lattice size L = 256 and 1024. Then in Fig. EJ^a), the fluctuation 
function F(t), defined in Eq. ((6]), is plotted for the depinning transition and its background, 
respectively. And basing on Eq. f|T6|) . the exponent 1/z — 0.749(5) and l/zj, = 0.677(3) are 
derived from the slopes of the dash lines. Here z is the dynamic exponent for the depinning 
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transition with A = 1.5, H c = 1.2933 and z& is for the background with A = 0, H c = 1.2933. 
Then the time evolution of the derivative d T In v(t, r) | T= o is displayed in Fig. [2(b) with open 
circles and its slope 0.729(4) is measured. In order to fit the numerical data well, a power 
law correlation is introduced, according to Eq. ( 1T7|) . With Eq. (fTUj) . 1/i/z — 0.735(3) can be 
measured. Following, in Fig. ??(a), the roughness function u) 2 (t) is plotted with open dots 
and their slopes are 1.701(6) for the depinning transition and 0.666(4) for the background. 
Then the pure roughness function Dcj 2 (t), defined in Eq. (fill) , is shown with pluses and the 
exponent 2(/z = 1.717(5) is obtained according to Eq. ( !T3|) . using this power law correction 
form again. In the same way, the exponent 2^/zi, = 0.649(4) is also derived. Here ( and 
denote the roughness exponent for the depinning transition and the background, respectively. 

Combing the exponents /3/vz, 1/vz, 1/z and 2(/z, one can calculate these critical expo- 
nents (3 = 0.295(2), v = 1.02(1), z = 1.33(1) and ( = 1.14(1) independent for the depinning 
transition. Then for its background, Zf, = 1.50(1) and = 0.487(5) are obtained. Fi- 
nally, in Fig. ??(b), the curves about the pure height correlation function DC(r,t) are 
plotted with open dots. For a sufficiently large scale r ^> £(;£), e.g., r = 256, the exponent 
2(/z = 1.701(7), derived with the scaling analysis in Eq. ( fl4l) . is well consistent with 1.717(5) 
obtained in Fig. ??(a). What's more, for a sufficiently small scale r -C £(£), e.g., r = 2, the 
curve also obeys a power law behavior and the exponent 2(( — Cioc)/z = 0.597(4) is also ob- 
tained. With z = 1.33(1) at hand, one can calculate that 2(i oc = (1.701—0.597) x 1.33 = 1.47. 
Then in the inset, the height correlation function DC(r, oo) vs. r is plotted for the depin- 
ning transition at the sufficiently larger time, e.g., t = 1024. And a special scaling form in 
Eq. (1151) is used to fit the data in the whole regime and 2Q oca = 1.46(3) is measured for the 
depinning transition, which is also well consistent with the exponent 2(i oc = 1.47 obtained 
before. 

Finally, the critical exponents characterizing the velocity of the domain interface, the 
roughness function and the height correlation function are extracted, and all the results 
are summarized in the last column of Table [B For comparison, some exponents, obtained 
by the driven QEW equation, the driven RFIM in the stead state, the discrete dislocation 
dynamics at the depinning transition and the forced-flow imbibition in columnar geometries, 
are also listed here. As was mentioned in the introduction, there are arguments in favor of the 
conjecture that the motion of a domain wall in a random-field system can be described by the 
driven QEW equation with quenched disorder [5^]. While our findings doesn't support this 
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conjecture, because there are at least 10% difference between these two models. Hence, they 
are in the different universality classes. While the theory values of these exponents for QEW 
have been derived by the function renormalization group(FRG) with z = 1.33, ( = 1.00 for 
the one-loop order and (3 = 0.31, = 0.98, £ = 1.43 for the two-loop order [291 ] . With 
this, one can see that some critical exponents in our work is well consistent with the theory 
value. Then these two measurements of the roughness exponent ( by the pure roughness 
function Du 2 (t) and the pure height correlation function DC(r,t) are in good agreement. 
And ( = 1.13(1) > 1 and (i oc = 0.743(6) < 1 mean that the interface growth process with 
intrinsic anomalous scaling and spatial multiscaling takes place. It is quite different with 
the interface growth process with superrough scaling and spatial single scaling in QEW with 
( = 1.23 > 1 and (i oc = 0.98 ~ 1 . The dynamic equation QEW is used to describe a one- 
dimensional driven elastic string in a two-dimensional disordered medium. Then the position 
of the interface can be defined uniquely and there are practically no overhangs, bubbles or 
droplets. So the single spatial scaling behavior is easy to understand. While overhangs is 
obvious in our work of RFIM, shown in Fig. [Jja). Hence, we guess it is overhangs cause 
this spatial multiscaling behavior, which is rarely in QEW. And additional simulation for 
the background in our work also supports this hypothesis. In our work, we confirm that the 
one-order transition takes place at H c = A, for the pinning- depinning process with A < 1 



26j. Hence, in our simulation of the background without disorder, a steady velocity state 



is obtained with the large driving field H ^> H c = 0. And a similar work has also been 
performed in the driven elastic string in a disordered medium [321 ] . To our surprise, the 
difference of the roughness exponent ( between these two models are neglectable within the 
errors, different form that at the depinning transition. And the overhangs vanish in the 
background with the finite velocity, which can be observed in Fig. []Ja). Hence, we can 
believe that overhangs lead to the main difference between this spin model and QEW. And 
( = 0.5 is expected, for the quenched disorder acts effectively as a thermal noise at the 



largest scale with (therm = 1/2 [34], 



35, 



611]- 



Then in the driven RFIM, the transition point H c and critical exponents /3, u, z, ( can also 
be measured in the steady state with finite-size scaling and finite-temperature scaling forms. 
Based on the short-time dynamics scaling form, it is more convenient and precise for the 
determination of both dynamic and static exponents as well as transition point. Since one 
does not suffer from critical slowing down. For instance, the transition point H c is reported 
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as 1.290(5) in Ref. [28[. While this curve with H = 1.290 visibly departure from the power 
law behavior, shown in Fig. QJb) . It means that the transition point H c is distinctly above 
1.290. What's more, in this work, we can investigate the roughness behavior carefully and 
systematically at the depinning transition. While for the RFIM in the steady state, it is 
too hard to measure the roughness exponent ( and dynamic exponent z. And in TableJH 
these exponents with + are obtained in the critical avalanches process, instead of depinning 
transition. Then the critical exponents of the depinning transition is also obtained for the 
discrete dislocation dynamics, which are quite agree with the result of our work. 

Experimental investigation of domain wall motion in magnetic systems takes place at 
finite temperatures. With the thermal activation, the sharp depinning transition is clearly 
rounded. In the creep regime with H <C H c , T > 0, the interface velocity does not vanish and 

And the local roughness exponent 



34 



43 



the scaling behavior has been investigated 

has been measured experimentally, shown in Tab. [Ill It is a litter larger than 



= 2/3, 
351 ]. influence 



obtained in the equilibrium state with H = 0,T > 0. Then in Refs.[27|, [28|, 
of temperature on the depinning transition of driven interfaces has also been shown with 
H = H C ,T > and the roughness exponent C — 1.25 is obtained, the same as the one at 



the depinning transition 



35]. And in the Ref. 34], ( = 1.26(1) can also be obtained, even at 



the driving field H well below H c . While in these experiment, H is not so small. Especially 



in Ref. 54J, H = Yl.lkA/m is used, quite near the critical field H c « 40kA/m. And these 
experimental result of the local roughness exponent are close to (i oc = 0.735(5), measured 



in our work. Hence, a crossover takes place from the equi 
states, as the driving filed is increased ; 



ibrium state to the depinning 



341 ] . And in Ref. 6l|], a similar crossover has been 



investigated as the temperature is decreased. 



IV. CONCLUSION 



In summary, with Monte Carlo simulation we have investigated the non-equilibrium crit- 
ical dynamics of the two-dimensional driven random-field Ising model for the depining tran- 
sition and its background at zero temperature. The dynamic scaling behavior is carefully 
analyzed, and a dynamic roughening process is observed. The transition point and critical 
exponents are totally listed in Table [H Based on the short-time dynamics, it is more conve- 
nient and precise to measure these exponents than the work in the steady state close to the 
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depinning transition. The conjecture that the motion of a domain wall in a random-field 
system can be described by the equation EW with quenched disorder is in doubt here. And 
intrinsic anomalous scaling and spatial multiscaling are observed in the interface growth 
process of our work, quite different from the superrough scaling and spatial single scaling in 
QEW. When the overhangs vanish with the large driving field, the difference between these 
two models disappears. Hence, we conjecture that overhangs lead to the main difference 
between QEW and RFIM at the depinning transition. What's more, the local roughness ob- 
tained in this work is supported by the results of experiment, shown in TabJITl But the effect 
of overhangs should be studied further. And it is also important to investigate the short- 
time dynamic behavior at lower temperature in RFIM and/or the transition between the 
relaxation, creep, slide, and switch regions with an oscillating driving field H = H exp(iut) . 
Finally, the techniques used in this paper can be also applied to similar dynamic systems, 
such as xy model. 
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Exponent QEW RFIM Dislocation This work 

v(t) H c 1.290(5)[28] 1.2933(2) 

P 0.33[33];0.33(2)[32] 0.35(4)[28];0.42(5)[62];0.31(8)[63] 0.30(5)[24] 0.295(2) 

v 1.33[33];1.29(5)[32] 1.00(5) [28] ;1.33[62j 1.05(5)f24] 1.02(1) 

z 1.5[33];1.54(5)[29] 1.32(4)f24] 1.33(1) 

u 2 (t) C 1.25[33];1.26(1)[32] 1-14(1) 

C(r,t) ( 1.2[64.];1.23(1)[45] 0.98(3)[24] 1.13(1) 

Cloc 0.92[64];0.98[31] 0.96(2)[24] 0.735(5) 

H->H C z b 1.5 [22] 1.50(1) 

( b 0.5[32J 0.5[62J 0.487(5) 

TABLE I: Summary of critical exponents are shown for the depinning transition (upper sector) 
and its background(lower sector), obtained with different techniques. Second column: result of 
numerical simulation obtained for the driven Edwards- Wilkinson equation with quenched disorder 
(QEW). Third column: result of numerical simulation obtained for the driven random field Ising 
model (RFIM) close to the transition point in the steady state. Forth column: result of numerical 
simulation obtained for the discrete dislocation model at the depinning transition. Then our work 
is shown in the last column, and the result of numerical simulation is obtained for the random field 
Ising model at the depinning transition with the short-time critical dynamic. 
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loc 



This work Ref. [5Q 
0.735(5) 0.69(7) 



Ref. 



54] 



0.78(1) 



Ref. 
0.7(1 



Ref. 
0.6(1) 



TABLE II: The local roughness exponent for the depinning transition in our work is shown here. 
In comparison, the result of experimental investigation of domain wall motion in magnetic systems 
are also listed in its creep regime. Refs. 



50h52|] work in the ultrathin Pt/Co/Pt films and Ref. 



works in Co2sPt72 alloy films. And an oscillating driving field H(t) is used in Ref. 
a dc driving field. 



54] 
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FIG. 1: (a) The time evolution of the interface velocity v(t) is displayed for different driving field 
H when the disorder is fixed as A = 1.5. For clarity, the curve with depinning field H c = 1.2933 
is shifted down by a factor 0.5. Dashed line shows the power law fit. Open squares corresponds to 
the different lattice size with L = 1024. (b) The function F(t) is displayed with open dots on 
a log- log scale for the depinning transition (above) and its background (below). And dashed lines 
show the power law fits. 
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t t 

(a) (b) 

FIG. 2: (a) The roughness function ui 2 (t) is plotted on a double-log scale with open dots. And 
pluses represent the pure roughness function Dui 2 (t). Solid lines represent power law fits with 
correction, basing on Eq. 1171 (b) The pure height correlation function DC(r,t) is plotted with 
open circles on a log-log scale, for r = 2,4,8, 16,32,64 and 256 (from below). Dashed lines show 
the power law fits , and solid line represents a power law fit with correction. In the inset, DC(r, oo) 
vs. r is plotted for the depinning transition at its maximum time t = 1024. Data are fitted with 
solid lines, according to Eq. [T5l 
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